Conductance of Tomonaga-Luttinger liquid wires and junctions 
with resistances 
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Abstract. - We study the effect that resistive regions have on the conductance of a quantum wire 
with interacting electrons which is connected to Fermi liquid leads. Using the bosonization for- 
malism and a Rayleigh dissipation function to model the power dissipation, we use both scattering 
theory and Green's function techniques to derive the DC conductance. The resistive regions are 
generally found to lead to incoherent transport. For a single wire, we find that the resistance adds 
in series to the contact resistance of h/e 2 for spinless electrons, and the total resistance is inde- 
pendent of the Luttinger parameter Kw of the wire. We numerically solve the bosonic equations 
to illustrate what happens when a charge density pulse is incident on the wire; the results depend 
on the parameters of the resistive and interacting regions in interesting ways. For a junction of 
Tomonaga-Luttinger liquid wires, we use a dissipationless current splitting matrix to model the 
junction. For a junction of three wires connected to Fermi liquid leads, there are two families of 
such matrices; we find that the conductance matrix generally depends on Kw for one family but 
is independent of Kw for the other family, regardless of the resistances present in the system. 



Introduction.- It is well known that for non- interacting 
electrons, the conductance of a narrow ballistic quantum 
wire is quantized in units of e 2 jh at low temperatures [112] , 
This remains true when electron-electron interactions are 
taken into account in the wire, provided that there are no 
sources of backscattering (such as impurities) and that the 
wire is connected to leads where there are no interactions 
[3H7]- Namely, if the wire is modeled as a Tomonaga- 
Luttinger liquid (TLL) and the interaction strength is 
given by the Luttinger parameter Kw, the conductance 
of a clean wire is independent of Kw- This breaks down 
if there are isolated impurities in a wire with interact- 
ing electrons; the impurity strengths then satisfy some 
renormalization group (RG) equations, and the conduc- 
tance depends on Kw and other parameters like the wire 
length, the distances between the impurities, and the tem- 
perature [8Ul0j . One can think of the impurities as giving 
rise to a resistance which leads to power dissipation, al- 
though this aspect is usually not highlighted in the liter- 
ature. There have been some studies of power dissipation 
on the edges of a quantum Hall system jllj and also at a 
junction of quantum wires due to the presence of bound 
states [12] . However, there has been relatively little discus- 



sion of the effects of an extended region of dissipation (a 
patch of resistance) within the framework of TLL theory 
or bosonization which is well-suited for studying the ef- 
fects of interactions between electrons [13] . Such a theory 
would have the benefit of combining the wealth of knowl- 
edge of TLLs with the classical notion of resistance. Fur- 
ther, a large amount of work has been done on junctions 
of several quantum wires theoretically [14fE5] and exper- 
imentally [24][25], and it would be useful to know what 
effect resistances in the wires have on the conductance 
matrix of such a system. A junction of three quantum 
wires with interacting spin-1/2 electrons has been studied 
in Ref. [55], and it has been found that some of the fixed 
points of the RG equations have different properties for 
the charge and spin sectors. In this context, we would 
like to mention the work in Ref. [27]. Here the effect of 
an extended region of inhomogeneity in a quantum wire 
has been studied, and it has been shown that this leads to 
weak backscattering which gives rise to a resistance which 
is linear in the temperature. Further, the resistances for 
the charge and spin sectors are different; the sum of the 
two gives the total resistance. 

In this paper, we will use the technique of bosoniza- 
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tion to study the effect of patches of resistance on the 
conductance of a quantum wire system with or without 
junctions. Our treatment will be classical in the sense 
that the resistance will be taken to be purely a source 
of Ohmic power dissipation; we will not consider the mi- 
croscopic origins of the resistance such as point impurities 
which can scatter the electrons quantum mechanically. As 
a result, the transport will be seen to be incoherent, with 
the resistance of different patches adding in series with no 
effects of interference; the incoherence also implies that 
RG equations will play no role in the analysis. Using the 
idea of a Rayleigh dissipation function [35J to model the 
resistance patches, we will obtain the equations of motion 
for the bosonic field whose space and time derivatives give 
the electron charge density and current respectively. We 
then use both a scattering solution of the equations of mo- 
tion [51|S] and a Green's function approach [3J to obtain 
the DC conductance G. Our analysis leads to several new 
results which are as follows. For a single wire, we cal- 
culate G when the Luttinger parameter K, the velocity 
v and the resistivity r all vary with the spatial coordi- 
nate x in the region of the quantum wire. The equations 
of motion enable us to numerically study the space-time 
evolution when a charge density wave of arbitrary shape is 
incident on the wire region. For the case of several quan- 
tum wires meeting at a junction, we model the junction 
using an orthogonal and dissipationless current splitting 
matrix M. For a three- wire junction, it is known that 
there are two families of M which have determinant ±1 
respectively When resistance patches are then introduced 
in each of the wires some distance away from the junction, 
the conductance matrix G of the system is generally found 
to depend on the matrix M as well as some of the param- 
eters mentioned above; this will be discussed in detail. 
For simplicity, we will restrict our analysis to the case of 
spinless electrons. 

Equation of motion.- We begin by studying a single 
wire with interacting spinless electrons. In the absence of 
backscattcring processes, the bosonic Lagrangian is given 

by 
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where K and v denote the Luttinger parameter and veloc- 
ity respectively; these parameters can vary with x within 
a finite region which we will take to be —L/2 < x < L/2. 
The Fermi liquid leads will be assumed to lie in the re- 
gions | a; | > L/2, where v = Vf and K = 1 arc constant; in 
the leads, the frequency and wave number of a plane wave 
are related as lj = vp\k\. The electron charge density n 
and current j are given in terms of the bosonic field as 
n = —edx^l^pK and j = edt<fi/ 'y/n, where e is the electron 
charge; these densities clearly satisfy the equation of con- 
tinuity dtJi + d x j = 0. The energy of the system is then 
given by 



We now introduce dissipation in the model through a 
Rayleigh dissipation function 



dx r j z 



(3) 



where the resistivity r can also vary with x but will be 
taken to be non-zero only within the region \ x\ < L/2. The 
function in eq. ([3]) contributes to the equation of motion 
as d/dt(5L/Sd t (t>)-5L/5<j>+5F/5d t (j) = [28], which gives 



4- dU - d x {^d x 4>) + -rd t cp = 0. (4) 



(Note that we have set h = 1, so that e 2 /(2ir) = e 2 /h). 
One can then show that the power dissipation is given by 



dE 
~dt 



dx d t 4> /8df 
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which is equal to — j 2 R in a steady state as desired; here 
R = f/^ dx r is the total resistance, and steady state 
means that j is independent of x (this follows from the 
equation of continuity and the fact that dtn = in a 
steady state). 

One can compute the conductance of the system in two 
ways. The first way is to consider the interacting and resis- 
tive regions as sources of scattering, as has been done for 
an interacting region in Refs. [5j[6]. We allow a plane wave 
with frequency uj to be incident on this region from the 
left, and compute the reflection and transmission ampli- 
tudes as functions of uj. The latter amplitude is related, in 
the limit uj — > 0, to the dc conductance a c i c . The second 
way is to compute the Fourier transform of the Green's 
function in imaginary time and hence the nonlocal ac con- 
ductance; this again gives o~dc m the limit u> — > [3J- We 
will use both these methods, the scattering method for a 
single wire and the Green's function for a junction of three 
wires. 

Transmission through a dissipative region. - To illustrate 
the scattering method for computing o^c [SJE], let us first 
consider a non-interacting system in which K = 1 and 
v = vp are independent of x, while r(x) = r$ for —a < 
x < a and elsewhere. This describes a dissipative region 
(—a, a) connected to leads on the two sides. For a plane 
wave incident from the left with k = lj/vf, the spatial 
part of the solution (j>k(x,t) = fk{x)e~ lujt is given by 



fk = e lkx + s k e~ lkx 
= t' k ,e ik ' x + s' k ,e 
= t k e lkx 



for x < —a, (6a) 
- lk ' x for -a<x<a, (6b) 
for a < x. (6c) 
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Using cq. ([3]), and the continuity of f k and d x f k at x = ±a, 
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we find that Sk and tk are given by 

1^2 _ ]Wgi2j7fca _ e -i2r)ka\ 

(1 + T]) 2 - (1 - rj) 2 e iAr >ka ' 
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where r\ 
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We now consider what happens when a ^-function 
charge density pulse is incident on the dissipative region 
from the left lead. At t = 0, the pulse is at x < —a with 
velocity +vp. This pulse is described by: 



4>{x,t) 



dk ^(x — XQ—VFt) 



2tt 



k- 



(9) 



so that n = eS(x — xq — vpt) and j — evp5(x — xq — vf£) 
for x < —a and t < —(a + xq)/vf (before scattering into 
the resistive region). For x > a, the corresponding current 
is given by 



evF 



dk 
2^ 



tk e 



ik(x — Xq— Vpt) 



(10) 



The nonlocal ac conductivity is given by a(x,Xo,t) = 
ej(x, i)/(27r), where x > a and xq < —a- The DC conduc- 
tance is given by the zero frequency limit of the Fourier 
transform, 



Ode 
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where x = (1 ~ 7 ?) 2 /(l + v) 2 - Using lim^o rj = \fCJk and 
limfe-s-o X — 1 — 4y / X7C, we find that 
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where R = 2ar$ is the total resistance. This expression 
shows that R adds in series to the contact resistance of 
27r/e 2 ; this is the expected property of resistance in a 
phase incoherent system. We obtain the same result for 
adc through the Green's function method outlined below 
for a three-wire system. 

The expression in cq. (|12[) can be derived for a general 
resistance profile r(x) as follows. For uj = 0, a solution 
of eq. is <f> = c. where c is a constant. Let us now 
look for a solution which is valid upto first order in uj and 
k = uj/vp, and which reduces to = c in the limit ui — > 0. 
Assuming that (f>k(x,t) = fk{x)e~ lu)t , where has the 
forms given by eq. (|6aj) and eq. (|6c| for a; < —a and > a 
respectively, we must have 1 + Sfc = tk = c to zero-th order 



in w as cj — > 0. Next, on ignoring the term of order w 2 in 
eq. ((!]), we obtain 



iuj— r fk = 0. 

7T 



(13) 



Replacing fk by the constant c = in the second term 
in eq. (fl~3|) (we can do this since that term has a factor 
of uj), and integrating that equation from x = —a — e to 
a + e, we obtain ikvp[tk — (1 — Sfe)] = —iLutk(e 2 R/ir), where 
i? = f ° q dx r(x) is the total resistance. Combining this 
with 1 + Sfe = ife , we obtain 



tk- 



1 



1 



2tt 



(14) 



from which the result for Ode = (e 2 /2-K)tk-+o follows. 

It is interesting to compare the evolution of a charge 
density pulse incident on a dissipative region in a non- 
interacting system with the evolution of the same pulse 
in a non-dissipative but interacting system (Kw ^ 1)- 
[The latter case was studied in Refs. [5j[6]. It was shown 
there that a series of pulses emerges on both sides of the 
wire, such that eventually the integrated pulse on the left 
(i.e., the total reflection probability) is zero, while the in- 
tegrated pulse on the right (i.e., the total transmission 
probability) is unity] . We have time evolved eq. Q numer- 
ically; a von Neumann stability analysis was performed to 
ensure that the numerical errors remain small. 



time = o.o 



time = 2.0 




Fig. 1: Evolution of a charge density pulse incident on a dissi- 
pative region (shaded), with lo = 0.1 (the width of the pulse), 
a — 0.25 and ro = 3. 

Figures [T] and [5] show the density profiles at different 
times for the purely dissipative and purely interacting 
cases respectively; in the first case, we have chosen K = 1 
and v — 1 everywhere, while in the second case, wc have 
chosen K = 1 and v = 1 in the leads, but K = 0.6 and 
V = 1.6 in the interacting region. In fig. [TJ one sees only 
one reflected and one transmitted pulse. The width of the 
reflected pulse (= 4a) is equal to twice the length of the 
dissipative region, providing us with an insight into the 
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Fig. 2: Evolution of a charge density pulse incident on an inter- 
acting region (shaded), with lo = 0.2 (the width of the pulse), 
= 1, Kw — 0.6 and vw = 1.6. 



nature of dissipation. Namely, a pulse gets reflected from 
each point in a dissipative region; this makes the width of 
the reflected pulse (i.e., the distance between waves being 
reflected from the left and right ends of that region) equal 
to 4a. On the other hand, when a pulse approaches an in- 
teracting region in which K is piecewise constant, it gets 
reflected only from the points of inhomogeneity, i.e., where 
dK/dx is not zero. Fig. [5] shows a series of reflected and 
transmitted pulses in agreement with the results obtained 
analytically in Rcfs. [5115]. 

Green's function calculation for three-wire junction. - We 
now consider a junction of three dissipative TLL wires as 
shown in fig. [3J (a), each of which contains three regions: 

(i) < Xi < Ln 7^ — the region around the junction 
where K(xi) = Kw, elsewhere K{xi) = 1, 

(ii) Ln < Xi < Li2 — a dissipative region where r(xi) = 
rio] elsewhere r(xi) = 0, and 

(iii) Xi > Li2 — semi-infinite leads. 

Here i labels the wires, and on wire i, the coordinate 
Xi runs from to oo, with Xi = corresponding to the 
junction point. The regions Xj > Lj2 model the two- or 
three-dimensional leads which are assumed to be Fermi 
liquids with no interactions between the electrons; hence 
we set K = 1 in those regions. 

Following Ref. [3], we write 
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in the linear response regime, where E u (xj) is the Fourier 
component of the electric field E{x'pt) on wire j, and 
o~ij,u(xi,%j) is the nonlocal conductance matrix. We then 
obtain 



i \ n i i \ 

\XijXj) yij.uj\Xi, Xj ), 



(16) 
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Fig. 3: (a) A schematic diagram of a three-wire junction with 
an interacting region close to the junction (red solid region), 
a dissipative region (blue dotted region further from the junc- 
tion), and Fermi liquid leads (brown solid region furthest from 
the junction), (b) The matrix elements Gij, with ij = 11, 12 
and 13, are shown as functions of Kw for the case of Mi with 
6 = 271-/3 and r i0 = for i = 1, 2, 3. 



where lu = —iu, and 



&ij,uj {Xi , Xj ) 



dr 

— {T*(j)i{xi 



-)&K-,o)> t 



(17) 

is the propagator of the bosonic field in imaginary time, 
t = it. The Green's function satisfies the equation 



v(Xi) 



K(xi) dx * 



2 — 



v{xi)K(xi) 



r(xi) 



(18) 



with the following boundary conditions: 
(i) Qij,a {%i j x 'j ) is continuous at Xi = x' 

a 



I/j2) and 



3 



Gij,uj (Xi, Xj)\ J, _ e 



j (where < x'j < 



v(xi) 



d Xi Gij,a{xi,x'A) are continuous 



dxi/viv r Qr Q < 



A.. ( fi 1 Xi/vw 

■f Ln(\ — dij) (vw is the velocity in 



(ii) Gij,a{xi,x'j) and 
at Xi = Lu and Li2, 

(iii) if Gij, a (xi, x'j) = 
Xi < mm.{x'pL il ) % 
the wire region), then B = —MA, where M is the current 
splitting matrix at the junction [T4H23] . 

The boundary condition in (iii) arises from the fact 
that the incoming and outgoing currents (and hence the 
bosonic fields) at the junction are related by the matrix 
M. Various constraints at the junction such as current 
conservation and unitarity of the evolution of the system 
in real time (i.e., no power is dissipated exactly at the 
junction) imply that each row and column of M must add 
up to unity and that M must be orthogonal. The possible 
M matrices are restricted to two classes parameterized by 
a single parameter 9 pP£M23) : (a) det(Mi) = 1 and (b) 
det(M2) = — 1, which can be expressed as: 



Mi 




and Mo 




,(19) 



where a = (1 + 2cos6>)/3 and 6(c) = (1 — cos# + 
(-)\/3sin0)/3. We note that (M 2 ) 2 = 1 for any value 
of 6; this relation will be used below. 
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Note that by introducing the orthogonal matrix M, we 
have made the simplifying assumption that there is no 
dissipation exactly at the junction. It is calculationally 
simpler to separate the junction, which governs how the 
incoming currents are distributed amongst the different 
wires, from the regions of dissipation which lie away from 
the junction. 

Solving eq. (|18p with the above boundary conditions 
and finally taking the limit ui — > + e, we get the following 
expression for the dc conductance matrix: 

G = - [t + M + K w (t~M)(t+-R)]- 1 

TT TT 

X [1-M], (20) 

where IR is a 3 x 3 diagonal matrix with Ru = R4 = 
fio(Li2 — Lii); note that Ri is simply the total resis- 
tance in wire i. The conductance matrix relates the out- 
going current Ii to the potential Vi applied in lead i as 
Ij = J2jGijVj. One can show in general that each row 
and column of G must add up to zero; the columns adding 
up to zero is a consequence of current conservation Ii 
must be zero), while the rows must add up to zero because 
each of the Ii must vanish if the Vj 's have the same values 
in all the wires. 

One can show that the conductance of a single wire 
given in eq. [12] follows from eq. (|20j) if we choose M = M2 
with 9 — 0, 27r/3 or 47r/3. For instance, if 9 = 0, wire 3 
decouples from wires 1 and 2 (so that Gij — if either 

1 or j = 3), while the conductance across wires 1 and 

2 becomes independent of Kw and is given by eq. (fT2"|) . 
with R = r 1Q (L 12 — + r 2 o(^22 — L21) being the total 
resistance in wires 1 and 2. 

Eq. (|20|) can also be derived in general using the equa- 
tion of motion approach in the uj — > limit in the same 
way as described above for the single wire case. We find 
that the precise profiles of K(xi), v(xi) and r(x;) in the 
different wires are not important; all that matters is that 
the values of K and v are given by Kw, vw a s Xi — > + e 
and by 1, vp as Xi — > 00, and that the diagonal elements 
of IR are given by Ri = J dxir{xi). 

Conductance for the M\ class.- In the M\ class, the 
case 9 = is trivial because Mi(0) = 1 and G = 0. Let 
us now consider other values of 9. We find that in general 
G depends on Kw, 9, and the resistances Ri = rio(Li 2 — 
Lii). [An exception arises for the case 9 = tt where we 
find that G is independent of Kw and depends only on 
the Ri. This occurs whenever M 2 = t which is true for 
M\{tt) and also for the M 2 class for any 9 as discussed 
below.] The dependence of G on Kw for the M\ class 
is to be contrasted to the case of a single wire where the 
conductance is independent of Kw [SHI] • In fig- El (b) , we 
show the matrix elements Gn, G12 and G13 as functions 
of Kw for the case of M\ with 9 = 2tt/3 and IR = 0. 

In the limit that Ri — > 00 (which is physically relevant 
when Ri >> 7r/e 2 ), we find that the conductance matrix 



takes the simple form 



R1R2 + R2R3 + R3R1 

(— R 2 — R3 R3 R2 \ 

R 3 -Ri -R 3 Ri , (21) 

R2 Ri —Ri — R2 / 

which is independent of both Kw and 9. Interestingly, the 
form in eq. (|2~T|) is exactly the same as that obtained for a 
classical system in which three wires with resistances Ri 
meet at a junction. If a potential Vi is applied to wire i, 
and the potential at the junction is Vq, then the outgoing 
currents are given by Ii = (Vq — V.)/Ri for all i. Using 
Kirchoff's circuit laws to eliminate Vq, we find that the 
conductance matrix relating Ii to Vj is given by eq. (|21j) . 

Conductance for the M2 class. - In this case the property 
M 2 = 1 combined with eq. (|2U| can be used to prove that 
G is independent of Kw for any choice of 9 and R4. The 
exact expression for G turns out to be 

e 2 3(1 - M 2 ) 
G - ~7 D ' 
where D = 2(qi + g 2 + Q3) + cos9(gi + g 2 - 2g 3 ) 

-V3 sin 0(01 - q 2 ), (22) 

where Qi = 1 + (e 2 /Tr)Ri. We can see that G does not 
depend on Kw- 

Time reversal invariance.- It is interesting to look at 
our results from the point of view of time reversal (T) 
invariance. There are two sources of T breaking in our 
system: 

(i) the presence of resistances (i.e., dissipation) clearly vio- 
lates T . This is evident from eq. (UJ which is not invariant 
under t — > —t. 

(ii) the current splitting matrix M is T invariant only if it 
is symmetric. This is because the outgoing and incoming 
currents near the junction satisfy I ou t = MIi n , while T 
interchanges Im and lout- We then see that ij„ = MI ou t 
is satisfied only if A/" 1 = M T = M. 

We observe that M 2 is symmetric and therefore T in- 
variant for all values of 9, while Mi is T invariant only if 
9 = or tt. A junction described by Mi can exist only 
if T is broken, for instance, by applying a magnetic field 
through the junction, assuming that this has a finite cross- 
section. 

The two sources of T breaking mentioned above are not 
related to each other since one acts at the resistances and 
the other acts only at the junction. However, we showed 
above that if M is symmetric, i.e., is T invariant, then 
the conductance G is independent of Kw, regardless of 
the values of the resistances. Thus there is a remarkable 
connection between T breaking at the junction and the 
dependence of G on the interaction parameter Kw- 

Discussion.- To summarize, we have presented a phe- 
nomenological formalism which allows us to study the ef- 
fect of resistive regions in a quantum wire using the lan- 
guage of bosonization. This enables us to calculate the 
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conductance of systems in which both the interaction pa- 
rameter K and the resistivity r vary with x. The bosonic 
equation of motion makes it possible to visualize what 
happens when a charge density pulse is incident on the 
resistive and interacting regions. Finally, by introducing a 
current splitting matrix M to describe a junction, we have 
extended the analysis to a three- wire system. We find that 
G depends on Kw for the class M\ (except for the spe- 
cial case with 9 = tt), but not for the class M2 in which 
case we have found an analytical expression for G. Thus 
we have generalized the well-known results of Safi-Schulz 
and Maslov-Stone to include systems with junctions and 
resistances. It may be possible to test our results experi- 
mentally by, for instance, varying Kyy and 9 by applying a 
gate voltage and a magnetic field near a junction of quan- 
tum wires [T8J|S0], and measuring how this changes the 
conductance matrix. 
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